Decay of quasiparticles in quantum spin liquids 
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Magnetic excitations are studied in gapped quantum spin systems, for which spontaneous two- 
magnon decays are allowed by symmetry. Interaction between one- and two-particle states acquires 
nonanalytic frequency and momentum dependence near the boundary of two-magnon continuum. 
This leads to a termination point of single-particle branch for one- dimensional systems and to strong 
suppression of quasiparticle weight in two dimensions. The momentum dependence of the decay rate 
is calculated in arbitrary dimensions and effect of external magnetic field is discussed. 
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A number of magnetic materials with quantum dis- 
ordered, or spin-liquid, ground states have been discov- 
ered in the past two decades. The well-known exam- 
ples include the spin-Peierls compound CuGeOsf 1 dimcr 
systems Cs3Cr 2 Brg£ and TlCuCl3^ integer-spin antifer- 
romagnetic chains^ and many others. Common prop- 
erty of all these systems is presence of a spin gap in the 
excitation spectrum, which separates a singlet, S = 
ground state from 5 = 1 quasiparticles. The low- 
energy triplet of magnons can be split due to intrinsic 
anisotropics or under applied magnetic field. Recent neu- 
tron experiments^ on two organic spin-gap materials 
PHCC and IPA-CUCI3 have revealed drastic transforma- 
tion of triplet quasiparticles undergoing at high energies. 
Upon entering two-magnon continuum, see Fig. ^ spon- 
taneous (T = 0) decay of a magnon into a pair of quasi- 
particles becomes possible, which leads to a rapid de- 
crease in the quasiparticle life-time in the former case^ 
and complete disappearance of the single-particle branch 
in the latter system^ 

Quasiparticle instability is well documented for an- 
other type of quantum liquid — superfluid 4 He. Pre- 
dicted by Pitaevskii nearly fifty years agOf£&& this 
instability was later confirmed by neutron scattering 
measurements. 10 ' 11 In liquid helium, interaction between 
one- and two-particle states is enhanced in the vicinity of 
a decay threshold by a large density of roton states and 
leads to avoided crossing: the single-particle branch flat- 
tens at energies below twice the roton energy and ceases 
to exist completely above that energy scale. 

The aim of present article is to investigate similar ef- 
fects in quantum spin liquids. We consider an isotropic 
spin system with a quantum disordered ground state, 
which is separated by a finite gap A from low-energy spin- 
1 excitations. A bare dispersion of propagating triplets 
is given by e(p). Bosonic operators t pa and ij, Q destroy 
and create a quasiparticle with momentum p in one of 
the three polarizations a — x,y,z. Interaction between 
one- and two-particle states is, generally, described by 
two types of cubic vertices shown in Fig. E^,b. In super- 
fluid 4 He presence of such particle non-conserving pro- 
cesses is determined by a Bose condensate, which ab- 
sorbs or emits an extra particle&i^ In quantum spin sys- 
tems with singlet ground states one finds a more diverse 



situation. Symmetry of a majority of low-dimensional 
dimer systems, like, for example, a dimerized chain^ an 
asymmetric ladder and a planar array of dimerSj 14 allows 
presence of cubic vertices. These should be contrasted 
with a symmetric ladder^ where one- and two-magnon 
states belong to sectors with different parity under per- 
mutation of two legs and, consequently, do not interact. 
There is also a suggestion that cubic vertices exists in 
a Heisenberg antiferromagnetic spin-1 chain beyond the 
nonlinear sigma-model description. 16 In the following, we 
shall assume presence of interaction between one- and 
two-particles states and consider its consequences for the 
dynamic properties of a spin liquid. 

In Heisenberg magnets, rotational symmetry in the 
spin space fixes uniquely the tensor structure of cubic 
vertices. Specifically, the decay vertex, Fig. has the 
following form: 

' k,q 

Conservation of the total spin during decay process re- 
quires that two created spin-1 quasiparticles must form 




FIG. 1: Schematic structure of the energy spectrum of a quan- 
tum magnet with spin-liquid ground state. Thin solid line is 
a bare spectrum e(p), full solid line is a renormalized single- 
particle branch, shaded region is two-particle continuum. The 
decay instability threshold is denoted by p c . 
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FIG. 2: Two cubic vertices of the decay (a) and the source 
(b) types and their contributions to the normal (c) and (d) 
and anomalous (e) self-energies. 



an S — 1 state, which is imposed by the antisymmetric 
tensor e a/31 . The vertex is, consequently, antisymmetric 
under permutation of two momenta r(k, q) = — T(q, k). 
Similar consideration applies to the source-type vertex, 
Fig. [2p, which is antisymmetric under permutation of 
any two of three outgoing lines. 

Important kinematic property of the energy spectrum 
e(p) is a type of instability at a decay threshold momen- 
tum p c , beyond which the energy conservation 



e(p) = e(q) +e(P - q) 



(2) 



is satisfied and two-magnon decays become possible. Ex- 
tremum condition imposed on the right hand side of 
Eq. (0) yields that two quasiparticles at the bottom of 
continuum always have equal velocities vq = v p _q, 
where Vk = Vk£(k). Then, the simplest possibility is 
that both momenta are also equal with Q = |(p + G), 
where G is a reciprocal lattice vector. We have ver- 
ified the above assertion for several model dispersion 
curve9 i7 i i8 ag we n ag £ or a f ew experimental fitsASiiS The 

common feature of all analyzed cases is presence of two 
regimes: for a small total momentum p of a magnon pair 
the minimum energy corresponds to equal momenta of 
two quasiparticles, whereas for p near the Brillouin zone 
boundary two particles of the lowest-energy pair have 
different momenta. Such a bifurcation is closely related 
to appearance of bound pairs of magnons at large mo- 
menta for one-dimensional spin systems^*2& Kinematic 
instability at the decay threshold in all considered cases 
corresponds, however, to a decay into a pair of quasi- 
particles with equal momenta. In the following we shall 
focus on such type of a model-independent instability. 

Let us consider the second-order contribution to the 
normal self-energy from the decay processes shown in 
Fig. Standard calculation yields 



d D q 



£u(c,p)= 2 y 



r(q,p-q) 2 



(2n) D lj - e(q) - e(p - q) + »0 ' 



(3) 



where D is the number of dimensions. We shall be in- 
terested in the behavior of energy spectrum at small 



Ap = p — p c and Alj = lu — e c and expand accord- 
ingly all functions under integral in Aq = q — Q, where 
Q = |(p + G). The main difference with Pitaevskii's 
analysis^ is that the antisymmetric vertex vanishes at 
the decay threshold such that 

T(Q + Aq, Q Aq) « 2Aq • V k r(k, q)| k , q=Q . (4) 

Including all dimensional constants in the vertex <?3 ~ 
|Vkr(k, q)| and performing angular integration for D > 
1 we obtain 



ft 

£ n (u;,p) ~ -g\ \ 
Jo 



q 



D+l 



dq 



V2 Ap — Au 



iO 



(5) 



where v 2 = Vk£(k)|k=Q and A is a lattice cut-off. The 
self-energy remains finite as Aw, |Ap| — > 0, though it 
contains an important nonanalytic contribution E(w,p). 
We shall absorb a regular part of Su(w, p) into e(p) and 
determine a new corrected spectrum from the Dyson's 
equation: G~ 1 (w, p) = lj — e(p) — £(u>, p) = 0. 

Two other second-order diagrams constructed from cu- 
bic vertices are presented in Fig. |2Ji,e. Correction to the 
normal self-energy shown in Fig. is an analytic func- 
tion of frequency and momentum and can be absorbed 
into e(p). The anomalous self-energy £i 2 (w,p), Fig. [2s, 
has the same dependence on lj and p as the one-loop dia- 
gram (|5J). Using the Belyaev's expression for the normal 
Green's function of a Bose gas&iS 



G 1 (lo,p) = u-£(p)-E n (u,p) 

Ei 2 (w, p)E 2 i(w,p) 



lj + e(-p) + £1.1 (-w, -p) 



(6) 



one can conclude that inclusion of the anomalous self- 
energy simply modifies a coefficient in front of the non- 
analytic contribution (JSJ without changing its functional 
dependence. Similar result applies to the vertex cor- 
rection from multiple two-particle scattering processes. 
Since a one-loop diagram is finite for lj = e c p = p c , the 
renormalized vertex §3 also remains finite and the vertex 
correction amounts to replacement g| — > (73173 in front 
of the integral in Eq. (JSJ. Thus, the following analysis 
based on analytic properties of the one-particle Green's 
function is not restricted to the lowest-order perturbation 
theory used in derivation of Eq. @ and is quite general. 

One- dimensional gapped spin systems: Let us begin 
with a stable region, where spontaneous decays are for- 
bidden: Au < 0, Ap > 0, such that {v 2 Ap - Aw) > 0. 
Separating the nonanalytic part in the integral JSJ and 
expanding e(p) ~ e c + v\Ap we obtain for the inverse 
Green's function 



G 1 (uj,p) = Aw — viAp — Ay v 2 Ap — Ac 



(7) 



with A = 7r<7§/2. The one-particle Green's function with 
such a dependence on to and p has been briefly discussed 



in Ref. 0. In the following we give more details relevant 
for the case of Eq. (7J . 
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Condition for a pole G 1 (u,p) — is transformed to 
quadratic equation, which yields 



e(p) 



i>i Ap - 



A 2 

2 



A 4 



+ A> 2 -i; 1 )Ap . (8) 



Away from the crossing point for Ap X 2 /(v2 — v%), 
slope of the single-particle branch coincides with the bare 
magnon velocity de(p)/dp « < 0, whereas for Ap — » 
slope changes sign to de(p)/dp = u 2 > 0. In the crossover 
region Ap ~ A 2 /(v 2 — V\) the quasiparticle weight Z, 



9G" 



9w 



= 1 



A 2 + v /A 4 + 4A 2 (i; 2 - t ; 1 )Ap 
4(^2 - v\)Ap 



(9) 

is continuously suppressed and vanishes, once the single- 
particle branch touches the two-magnon continuum. 

On the opposite side of the crossing point, Ap < 0, 
Aw > 0, and (Aid — V2Ap) > 0, the nonanalytic 
part of the self-energy becomes purely imaginary S = 
— iAi/Aw — V2Ap. Formally, condition G _1 (w,p) = 
yields after transformation to quadratic equation the 
same solution as in the stable region. For small nega- 
tive Ap the root (JSJ is real and satisfies Aid — v\Ap < 
and Aid — V2Ap < 0. The latter of the above two relations 
contradicts the assumption made in the beginning of the 
paragraph for the decay region, whereas the former is 
inconsistent with zero of G _1 (w,p) in the region, where 
v 2 Ap— Aw > 0, sec Eq. (JJJ. Thus, inside the continuum, 
a physical pole of the one-particle Green's function reap- 
pears only at a finite distance from the crossing point, 
where Re {Aw} > v 2 Ap or 



Ap < 



A 2 



2(«2-«l) 



(10) 



Disappearance of singlc-magnon excitations inside the 
continuum is quite similar to the termination point in 
the energy spectrum of superfluid ^eii&^iiLII In both 
cases the single-particle branch approaches tangentially 
the boundary of the two-particle continuum. 

Further away from the lower edge of the continuum in 
the region, where e(p) « e c + V\Ap, one can find for the 
inverse quasiparticle life-time: 



7p ~ 2A v / (^2-wi)|Ap| 



(11) 



Such a square-root dependence of the decay rate follows 
also from a simple golden rule consideration^ We have 
seen, however, that the full dependence of the self-energy 
on frequency and momentum must be kept in the vicinity 
of the decay threshold. Besides, if the cubic vertices have 
no special smallness, i.e., A is of the order of the band- 
width of e(p), the perturbative regime, where Eq. 
applies, is never realized. Such a situation apparently 
occurs in the quasi one-dimensional gapped spin system 
IPA-CUCI3, where no trace of single- magnon excitations 
has been observed inside the continuum. 6 



Two-dimensional gapped spin systems: Below decay 
threshold in the stable region, where (v 2 Ap — Aw) > 0, 
the nonanalytic part of the self-energy becomes 



£(w, p) = A(v 2 Ap - Aw) In 



R 



v 2 Ap — Aw 



(12) 



where A = <?f/2 and R is an energy cut-off. Though, 
a pole of the Green's function cannot be found analyti- 
cally, it is easy to verify that solution of G _1 (w, p) = 
exists all the way down to the crossing point. The quasi- 
particle branch touches tangentially the boundary of the 
continuum V p e(p) = v 2 for p = p c , whereas the quasi- 
particle weight is gradually diminished and vanishes at 
the boundary. 

Above the decay threshold [(v 2 Ap — Aw) < 0] the 
self-energy acquires imaginary part: 



E(w, p) = -A(Aw - v 2 Ap) ( In 



R 



\ Aw — v 2 Ap 

(13) 

In the perturbative regime, where e(p) » e c + viAp, the 
decay rate of magnons grows linearly inside the contin- 



7 P ~ 2ttA(vi - v 2 )Ap 



(14) 



In the close vicinity of the crossing point we write instead 
e(p) w v 2 Ap + a with Re a > and transform equation 
on the pole of the Green's function to 



a 



(vi - v 2 )Ap 
A 



b > 



(15) 



With logarithmic accuracy the solution of the above 
equation is 



Re a — 



Im a = 



-ixb 



In 2 R/b 



In R/b ■ 



(16) 



Single-particle branch of a two-dimensional gapped spin 
system can be, therefore, continued inside the decay re- 
gion. However, the quasiparticle weight is suppressed in 
a range of momenta near the crossing point: 



(vi - v 2 )Ap| ~ Re~ 1/X 



(17) 



Suppression of the quasiparticle peak has been exper- 
imentally observed in two-dimensional quantum disor- 
dered antiferromagnet PHCC£ 

Three-dimensional gapped spin systems: The nonan- 
alytic part of the self-energy is given in this case by 



S(w,p) = -A(v 2 Ap- Aw) 3 / 2 , 
S(w,p) = -iA(Aw-v 2 Ap) 3/2 , 



Aw < v 2 Ap , 
Aw > v 2 Ap , 



where A = irg^/2. Such a weakly nonanalytic term can be 
treated perturbatively near the crossing point. In partic- 
ular, the quasiparticle weight is not suppressed as p — > p c 
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and e(p) does not change its slope at the crossing. The 
magnon decay rate in three-dimensional case is 

7p ~ 2A[( Vl - v 2 )Ap] 3 / 2 . (18) 

An experimentally relevant question is how an external 
magnetic field affects the decay processes. If the Zeeman 
energy is smaller than the spin gap g^BH < A, the role of 
an applied field reduces to splitting a triplet of low-energy 
excitations. This can, in principle, lead to breaking the 
resonance condition (J2J). To study such a possibility one 
should transform from the vector basis t pa to states with 
a definite spin projection on the field direction: 

t p o = t P z , t p ± = {t P x T Upy) . (19) 

The decay vertex (QJ taken in the new basis has the fol- 
lowing form 

^AA^ = (iA- iA-K* (2°) 

4- (t f f f - ^ /t \ f , ( f t A _ t t f t \ f 

In all decay channels a destroyed one-particle and a 
created two-particle states experience the same energy 
shifts. Therefore, the decay threshold momentum p c 
does not change for all three splitted branches of the 
spin-1 excitations. 

Once the Zeeman energy exceeds the triplet gap, H > 
H c = A/g/iBi a Bose condensation of magnons takes 
place in D > li2i The ground state acquires a nonzero 
magnetization and a long-range order of transverse spin 



components. Such a canted antiferromagnetic structure 
opens an additional channel for spontaneous decays of 
low-energy excitations in the gapless branch, similar to 
the prediction made for ordered antiferromagnets in the 
vicinity of the saturation fields 

Intrinsic magnetic anisotropics, which become impor- 
tant in materials with spins S > 1, can also affect the de- 
cay processes. Anisotropy changes differently dispersion 
of triplet excitations and modifies the resonance condi- 
tion (J2J). The decay rate of spin-1 excitations depends, 
then, on a polarization, as was experimentally observed 
in bond-alternating quasi-one-dimensional antiferromag- 
net NTENPiS^ The tensor structure of the decay ver- 
tex Q can be also modified due to the absence of spin- 
rotational symmetry. Dynamic properties of spin liquids 
with different types of anisotropics deserve further theo- 
retical investigations. 

In conclusion, by studying analytic properties of the 
one-particle Green's function near the decay threshold we 
have found that crossing of a single-particle branch into 
two-magnon continuum is described solely by a growing 
line-width of magnons only for three-dimensional quan- 
tum spin liquids. In two dimensions there is, in addi- 
tion, strong suppression of a one-magnon peak near the 
crossing point, whereas in one dimension a single-magnon 
branch terminates at the continuum boundary. 
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